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Abstract 

We study the time dependent transitions of quantum forced harmonic oscillator (QFHO) in 
noncommutative R^'^ perturbatively to linear order in the noncommutativity 9. We show that 
the Poisson distribution gets modified, and that the vacuum state evolves into a "squeezed" state 
rather than a coherent state. The time evolutions of uncertainties in position and momentum in 
vacuum are also studied and imply interesting consequences for modeling nonlinear phenomena 
in quantum optics. 

1 Introduction 



There seems to be a growing consensus among physicists that our classical notion of spacetime has 
to be drastically revised in order to find a consistent formulation of quantum mechanics and grav- 
ity [IHS]. One possible generalization that has attracted much interest is that of noncommutative 
Moyal spacetime [H-Hl]. In situations where the time coordinate remains commutative, i.e., only the 
spatial coordinates do not commute with each other, the quantum theory is conceptually straightfor- 
ward (but nonetheless may display novel phenomena) |T5H2T]. In this article we will concentrate on 
■rj" ! understanding some implications of quantum mechanics with time-space noncommutativity, specifi- 
^ ically we will work with the Moyal plane Hg'^. We will use the formalism of unitary quantum mechanics 
,—1 on this space as developed by Balachandran et. al. [22] (see also [23]). 

^ ■ When time and space do not commute with each other it is not unreasonable to expect that the 
kjp .dynamics of the time dependent processes get altered. We will verify this explicitly in the context of 
!a simple model of the forced harmonic oscillator (FHO) with the forcing term switched on only for 
'a finite duration of time. In the commutative case this is a much studied model. We will compute 
deviations from the commutative case to leading order in 6. These deviations suggest that time- 
space noncommutativity can capture certain nonlinear effects seen in quantum optics. This article is 
organized as follows: In section [2] we will briefly review the formulation of unitary quantum mechanics 
on R^'^ [22] • In section [3] we will solve the problem of the FHO perturbatively in 6 and compute 
corrections to the transition probabilities between simple harmonic oscillator (SHO) states. These 
corrections suggest the noncoherent nature of the time-evolved vacuum state and are the reminiscent 
of those seen in nonlinear quantum optics [21]. To flesh out this analogy better we study the time- 
evolution of uncertainties in position and momentum in section HI Encouraged by these results we, 
in section [5l suggest a correspondence between the nonlinearity in quantum optics and the quantum 
mechanics on Rg'^. We conclude with a summary of our results in section 
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2 Unitary Quantum Mechanics on rJ'^ 

The noncommutative space R^' ^ is described by the coordinates f^'s satisfying 

[Xf,, Xu] = iOe^u with e^^ = -e^^ and Eqi = 1, (1) 

where n and u can take values 0,1. Without loss of generality we can take 6' > 0, as its sign can 
always be flipped by changing xi to —Xi. Let AeilO-'^) be the unital algebra generated by Xq and Xi. 
We associate to each a ^ Ae (R^'^), its left and right representations and a^: 

a^P = d/3 , a^P = Pa , /3 e Ae (R^'^) . (2) 

Unless stated, we work with the left representation. 

For a quantum theory, what we need are: (1) a suitable inner product on Ae (R^'^); (2) a 
Schrodinger constraint on AeilO-'^); and (3) a self-adjoint (with respect to the inner product de- 
fined) Hamiltonian H and observables which act on the constrained subspace of Ae (R^'^). 

1. The Inner Product: 

There are several suitable inner products and they are all equivalent to each other. One example is 
given here: We associate a symbol 0:5 corresponding to each a = J (fka{k)e^^^^^e^^^^^^ G Ae{B}'^) as 

as{xo,xi) = j d^ka{k)e"'''''e""""> . (3) 

Note that xq and xi and hence as are purely commutative. The inner product is defined as 

(^a,(3j ^ = j dxia*s{t,xi)(3s{t,xi) . (4) 

2. The Schrodinger Constraint and time evolution: 
The operators Pq and Pi, given by 

d 1 d 1 

2-— = Po = --adxi, -i-— = Pi = --adxo, (5) 
oxq 9 0x1 9 

generate time and space translations respectively. The Hamiltonian H, in general, may depend on 
Xi,Xq and Pi- The possible dependence of xf and Xq can be bypassed by 

xf = 9P0 + , x^ = 9Pi + x^. (6) 

Also, there is no dependence on Pq assumed in the line of the commutative case where there is 
never such dependence of H on idxQ for ^ = 0. Now note that the inner product (jl]) has an explicit 
dependence on the parameter t and hence there exist more than one null vectors with respect to this 
inner product (actually any vector which vanishes at Xq = t is a null vector). But this fact need not 
bother us as we are only interested in those states that satisfy the Schrodinger constraint 

Po - ^) = 0. (7) 

It is easy to see that now there are no non-trivial null vectors. The Hamiltonian H depends on xf , x^ 
and Pi. Since x^ commutes with xf and Pi we will choose x^ as "time". 

It is easy to write down the formal solution of the Schrodinger constraint and find the time 
evolution. The time evolution is given by Xq — )■ Xq + t (or equivalently by x^ — )■ x^ + r). Thus the 
amount of time-translation is always commutative, though the time-operator itself is noncommutative. 
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The time evolved wavefunctions satisfying the Schrodinger constraint are of the form ■iIj{xo,xi) = 
U (f^, T/) x{xi), where 

(8) 

3. The Spectral Map: 

Consider a time-independent Hamiltonian H = ^ + V{xi). The corresponding commutative Hamil- 
tonian is H = — ^(^i)' with eigenfunctions ipE (xq^Xi) = (pE{xi)e~^^^° and eigenvalues 

E. The spectrum of the corresponding noncommutative H will be given by "ipE = e~'^^^0(pE{xi) = 
ipE{xi)e~''^^° with the same eigenvalues E as HipE{xi) = EipE^xi). Here (Pe{xi) has been obtained 
by replacing xi with xi in {Pe{xi). 



U{x^,Ti) = (Texp 



drH 



TT 



3 QFHO in and their Transition Probabilities 

Let us recall the dynamics of a QFHO in ordinary spacetime. The Hamiltonian of this system is given 
by 

Hit) = + -m^uj^x^ + f{t)x + g{t)p, (9) 
IniQ 2 

where mo is the mass of the particle and u is the angular frequency of the oscillator. We are interested 
in real functions obeying 

/(t),^(t) = for t ^ ±oo. (10) 

At t — )• — oo the Hamiltonian is simple harmonic and we assume the system to be in one of the 
eigenstates of this SHO Hamiltonian. At t — )■ oo the Hamiltonian again becomes simple harmonic 
and we try to find the probability (the transition probability) for the system to be in any arbitrary 
eigenstate of the SHO Hamiltonian subjected to the fact that the system was in some already given 
eigenstate at t — )■ — oo. For this what we do is the following: 

• First we assume our system to be in an eigenstate ipnix) at t = tj — i- — oo. 

• The state (pnix) evolves under the SHO Hamiltonian from t = ti ^ —oo to t = Ti. 

• At t = Ti the interaction gets switched on. 

• The system then evolves under the full Hamiltonian (jH]) from t = Ti to t = T2. 

• At t = T2 the interaction gets switched off. 

• The system again evolves under the SHO Hamiltonian from t = T2 to t = tj ^ 00. 

• We find the inner product of the final state we get at t = — )■ 00 with the eigenstate (j)m{x)- This 
gives the Transition Amplitude Amn while its absolute square gives the Transition Probability 
P 

mn- 

The generalization of the above Hamiltonian in Rg'^ is 

p^ 1 1 ^ ^ 

H = + -mou'^xj + -[/(xo)xi + xif{xo)] + g{xo)pi = Hq + Hi, (11) 

with 

1 - 1 

Ho = + -mou'^xl, Hi = -[/(xo)xi + xif{xo)] + g{xQ)pi. (12) 



As xq and pi commute with each other, the ordering does not matter in the last term. 

To define the transitions for the above Hamiltonian consider the time evolution by an amount r. 
The functions /(xq) and g{xo) have the properties of vanishing in the far past and the far future, i.e., 

/, g{xo + r) — > as r — 7- ±00. (13) 

We shall find the transition probabilities {Pm,n) for an SHO state "n" at initial time (r — t- —00) to 
go to some other SHO state "m" at final time (r — )■ +00) after evolving under the Hamiltonian f llip . 
The Spectral Map tells us that the energy spectrum of the SHO Hamiltonian in R^'^ is same as that 
of the commutative one, i.e., 

E^ = hw(^n + ^y M^o, xi) = </>„(xi)e-'"("+^)"°, (14) 

where 0n(xi) is the eigenfunctions of the commutative SHO Hamiltonian. The orthonarmality (apart 
from a phase factor which comes because of the time evolution) of the eigenfunctions ipn{xo,Xi) with 
respect to the inner product defined in section [2] can easily be checked. The transition probabilities 
for our problem can be found by computing the same for the commutative Hamiltonian obtained after 
replacing 

e , , 

a;o-^-i:P + ^? xi x, pi p (15) 
h 

in the Hamiltonian ( ITTj) . Here t has come in place of the "time" x^ which commutes with x\ and 
To linear order in Q we obtain the following commutative Hamiltonian 

H{t) = Ho + Hj{t) = Ho + Hjoit) + OHnit), (16) 

with 

Ho = £-o + l^^ow^x^ = hw (ata + i) , Hjo = f{t)x + g{t)p = z*{t)a + z{t)a) 



Hn = i {-g'{t)p^ - \f'{t){xp + px)) = iyf^^ (^z*'{t)a^ - z'{t)a^^ + i^J^^g' {t){2a^a + I] 
The function z{t) is related to f{t) and g{t) as 



(17) 



z{t) = ^l^^{f{t) + zmoUjg{t)) 



Also, a and are the annihilation and creation operators respectively defined as 



at = ^ x-^^ p = ^J^{a^-a) 



The nonlinearity in the Hamiltonian flTB]) is purely due to the noncommutativity. This provokes us 
to model certain types of nonlinear phenomena in quantum optics by the noncommutativity between 
time and space coordinates. This analogy will be further studied in section [51 Let us now continue 
with calculating the transition amplitude which is given by 

Am,n{tf,T2;Ti,U) = {<PmH{tf,ti)Uo{tf,T2)U{T2,Ti)UoiT,,ti)\(j)n), (19) 

where Uo{t', t) and U{t', t) are the time evolution operators from time t to time t' for the Hamiltonians 
Hq and H(t) respectively, i.e., 

Uo{t',t) = e-^^»(*'-*), f/(t',t) = T \e-TJt' , (20) 
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the latter one being the time-ordered exponential. This gives 

Am,n{tf,T2;T,,U) = e^[^™(^^-*^)+^-(*-^i)l(0^|f/(T2,Ti)|0„). (21) 
The state \ip{t)) = U{t,Ti)\(f)n) evolves according to the Schrodinger equation for the Hamiltonian 

m 

th^-Ho] \m) = Hiit)\m), (22) 



dt 

with the initial condition \ip{t = Ti)) = If we define the Green's operator function G{t,to) as 

zh^-Ho^G{t,to) = S{t,to), (23) 
then solution of the Schrodinger equation fl^ will be 



m)) = m)) + / G{t,to)Hj{to)mto)), 

J — oo 

which in turn gives the Born series 

/ + 00 
dtoGit,to)Hi{to)mo)) 
-oo 

/+00 r+oo 
dto / dt,Git,to)Hi{to)G{to,h)Hjihmh)) + .... (24) 
■oo J — oo 

Here \4>(t)) is the solution of the homogeneous equation {ih-^ — Hq) |0(t)) = 0, which is nothing but 
the Schrodinger equation for SHO. G has been found in the Appendix |X] (see fl77|) ). Note that the 
G-function in the expression of the G restricts the integration over tj in f l2^ within the limit of — cxd 
to tj-i (t_i = t). Thus, at t = Ti the integrations are only in the intervals when the interaction 
was switched off, i.e., Hj = 0. Hence, we get \ip{t = Ti)) = |(^(t = Ti)) = |0„). The solution of the 
homogeneous part \4>{t)) with this initial condition is 

|0(t)) = e-i^"(*-^^)|0„). (25) 

Now, putting ([21D with t = T2 for [/(Tg, Ti)|0„) in we get 

00 

Ar^^n{tf,T2;T,,U) = Y,Bj{tf,T2;T,,ti), (26) 

j=0 



with 

'-'m.n 

Bj{tf,T2]Ti,ti) = dtQj_^ dti...J_^ dtj_iFl^^{tf,T2]to,ti, ...,tj_i;Ti,ti 

for j = 1, 2, .... Here 

Fi^^itf, T2; to, ti, Ti, t,) = {-lY e(T2 - to)e{to - U)...Q{t,_2 - h-i) 

eii''--E^)*^{ct)rn\Hr\to)Hr\ti)...Hr\t,.,)\(l)n). 



(27) 



(28) 



The _f/'*"*'s are defined as 

Hf^^it) = e^^o*iJ(...)(t)e-*^o* (29) 
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Separating the 6'-dependent and independent part we get 

A^,4tf,T2;T,,U) = ei(^"-^-)*^(0„|[AW(T2,Ti) + M«(T2,Ti)]|0„), 

with 

+... 



(30) 



(31) 



T 



+.... 



(32) 



The above expression for A^^\T2^Ti) can be simphfied to (see Appendix [B|) 

/i«(r2,ri) = --A(°)(T2,Ti) y f/to [AW(to,Ti)] (to)^(°H^o,ri). (33) 

A'^'^\t^t') with arbitrary arguments is defined in fl80p . Putting this in equation fl3UI) we get 



A straightforward use of the identity 

A 



■ poo 

I--eJ dto[A^''\to,T,)]'''H}fito)A^'\to,T,] 



(34) 



A3 



^AA^^-AA = B + ^[A,B] + '- [A, [A, B]] + [A, [A, [A, B]]] + ... (35) 

gives all the if'"*'s. Also, to get rid of the time ordered exponentials we follow the discussions given 
in pages 338-340 of [25]. This finally gives the expression of the transition amplitude as 



+/32* V^5:T1I^™,„+i(0 + P3^/n{n-l)D^^n-2{0 
+13; ^{n + l){n + 2)D^,^+2 (0 } 



with 



dre'^^^zir) 



— oo 
oo 



— / dn I dT2 [z*(ri)z(r2)e-*-(^^-^^) - ^(ri)^*(r2)e^"(^-^^)] = real 



(36) 

(37) 
(38) 



oo ^ — oo 

+ 00 



drg'irM^, 



— oo 

+ 00 



h 



i rriohjjj 



+00 



mou / drg'{r)C{r) + 



2i ImQhuj 



— oo 

+ 00 



h 



i mQhw 



+00 



dTz*'(T)e 



-2iujT 



dr [z*'{T)^\T)e-^'^^ - z\T)C\r)e^''^^] (39) 
rfrz*'(r)^(r)e^2i..r ^^q^ 

(41) 



Here the function C,{t) is given as 

m = -\f_ dre'-^ziT). (42) 
-Dm,n(0's are the matrix elements of the displacement operator D{^) = e~^*"+'^'^^ given by [26] 




(43) 



L,^(x) are the associated Laguerre polynomials. Also, the limits of the integrations have been extended 
to — oo and oo as the integrands are zero in the extended region. The transition probability is given 
by 

= |A™,n(t/,T2;Ti,ti)|2 (44) 

as usual. The arguments have been omitted as the transition probability does not depend on the 
times tf,T2;Ti,ti. 



3.1 n = 

For the initial state \4>o), the transition amplitude is 



A^,o(t/,T2;Ti,t,) = e^/^ei(^o-^™)*»e- 



1 - - PIC + PlCn + m^m* - 2/33T') 



+m(m-l)^/33T'} 



and the transition probability becomes (upto linear order in 



Pmfl — I ^m,o {if ,T2;Ti,ti 



2m 



ml 



1 + -6{Ai + mA2 + mim - l)A^} 

h 



with 



(45) 
(46) 

(47) 



Note that as m — >■ cx3, the ^-correction starts dominating and in this case the expansion upto linear 
order in 6 is no more meaningful. Hence, the above result is valid only for those m-values which are 
far smaller than l/\/moOj6 (in the unit ^ = 1). For 6* — )■ 0, the transition probability becomes the 
well known Poisson distribution as expected. 

As a specific example let us work with the functions f{t) and g{t) of the form (see Figure [1]) 



m 

9{t) 



/o[e(t + r)-0(t-r)] 
go [e{t + T)-e{t-T)] 



;T > 0. 



(45 



For these functions we get 



Ai 
A2 
As 



2/0 

mod;' 



(/q + niQU go) sin uT cos2LdT 



mougo sin 2uT — fo cos 2uT 
mlu'^go cot uT 
f^ + mlu'^g^ ■ 



(49) 
(50) 
(51) 
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Now the choices mo = 1, a; = 1, /o = Vb, 
the following Poisson distribution: P, 



m,0 



90 

-20 20 



a/S and T 



I in commutative case 



0) give 



modifies to P, 



m.O 



-20 20'^ 



e while for nonzero 9 the probability distribution 

[1 + 2Q^^{m — 20)]. The ^-correction becomes of the order of the Q- 

independent part when m approaches the value miff) = ^20 + ^^7^^^ Hence, our result is valid only 

in the region where m < fh{9). Note that 171(9) ~ | rather than ^ because the A3 is identically zero 
for the choices taken. We choose 6 = 0.01 (m{6 = 0.01) ~ 42) and get 



m,0 



-20 



20' 



ml 



■[l + 0.02\/5(m-20)]. 



(52) 



This deformed distribution along with the Poisson distribution is shown in Figure |2j Such deformation 
of the Poisson distribution suggests that the vacuum does not evolve to be a coherent state anymore. 
To explore this further let us look at the time-evolution of position and momentum uncertainties. 



(1) 




^ fit), git) 








fo, go 










T 


+ 


T t 



0.08 



0.06 



0.04 



0.02 



0.00 



• — • — • Poisson 
Distribution 



Modified 
Distribution 




Figure 1: The behaviour of functions fit) and git) with t 

Figure 2: The modified distribution {9 = 0.01) along with the Poisson distribution {6 = 0) for itlq = 1, 
u = 1, fo = V5, go = V5 and T = | 



4 The time evolution of A^; and 

The expectation value of any operator (9 in a state ipixQ,Xi) at any time t is defined to be 

{d)t=[^,d^)^. (53) 

Also, 

{d)t+r= (^,0^)^^^= (v^(xo + r,xi),C'7/;(xo + r,xi))^. 

Thus the time evolution of the expectation value of an operator is given by that of the state in which 
it is being calculated. For the QFHO in R^'^ the time evolution of any operator O will be given by 

) , (54) 



H{t),d 
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where H{t) is the Hamiltonian f|T6|) . The uncertainty in any observable O is given by 



(55) 



Thus the evolution of and is 



dt^P 



-2mou^ (lixp+px) - {x){p)) + frmi 



Defining 



Axp = -{xp + px) - {x){p), 



(56) 



(57) 



we find the following first order coupled equations 



d_^2 

dt^x 

dt^P 
A A 



-2mou;2^.p + f /'(t)A2 



(5^ 



As the initial state is the vacuum, the initial conditions for the above are 



A2(t ^ -oo) 



h 



2moUJ 



Al{t -oo) = A,p{t -oo) 



0. 



(59) 



Our strategy for solving these equations is simple. We do so perturbatively in ^. A straightforward 
computation gives 



A..(t) 
Apit) = 



h 9_ I mp 

2m(,u] 2 V 2fia) 



2 2 V 2fi 



+ dTsm{2u{r-t)}f{T) + 2u' f dr cos{2uj{r - t)}g{r 

fit) + tJ dTs\Yi{2u{T - t)}/(r) + 2uj^ j dr cos{2u{t - t)}g{T 

J — oo J —oo 



2 

mo 



g{t) + 2u dr cos{2u{t — t)}f{T) — tuqu'^ dr sm{2uj{T — t)}g{T 



A,p{t) = 9 



The fundamental uncertainty product (to linear order in 6) is 

A,(t).Ap(t) = J. 



(60) 



(61) 



Thus the vacuum state evolves to a "squeezed state" rather than a coherent state as in the commutative 
case [27]. The uncertainties in the commutative case depend only on the product mow. But, their 
^-corrections change with u even if niQU is kept constant. Also, the squeezing effect is oscillatory in 
time as is obvious from the ^-dependent terms in (1601) . For the specific forms of f{t) and g{t) of 
we get 



A,{t) 



h 



2mQU} 



:t<-T 



h 



2mQU 



\— [2/o cos{2u;(t + T)} + mo^^ow sin{2a;(t + T)}] ; -T < t < T 



h 



2mou] 



2hmouj 

2 V 2H^ [2/o {cos{2u{t + T)} - cos{2a;(t - T)}) 

+ mogou (sin{2u;(t + T)} - sin{2u;(t - T)})] 



(62) 



■t>T, 



9 



Ap(t) 



and 



mphui 
2 

+ I [2/o cos{2u;(t + T)} + mo(?oW sm{2u{t + T)}] 

+ f [2/o (cos{2a;(t + T)} - cos{2a;(t - T)}) 

+ moc/oc^ (sin{2w(t + T)} - sin{2w(t - T)})] 



I [2/o sin{2u;(t + T)} - motugo cos{2u;(t + T)}] 



I [2/o (sin{2a;(t + T)} - sin{2a;(t - T)}) 



—rrioLugo (cos{2a;(t + T)} — cos{2a;(t — T)})] 
For the same choice of parameters as in the last section we get 



]t<-T 
■-T <t<T 

■t>T 

]t<-T 
~T <t<T 

\t>T. 



(63) 



(64) 



1 

V2 



■t < -- 



Ax(t) 



^ + 0.01 J| (cos 2t+ i sin 2t) ; -f < t < f 



(65) 



^ + 0.01^/1 sin2t ;t>f 



2 ' 



1 

71 -0-01 



(cos 2t + i sin 2t) 
I sin 2t ; t > f 



(66) 



(67) 



and A^p(t) = <{ 0.01^5 (| cos 2t - sin 2t) 
-0.02\/5sin2t ;t>f. 

Figures [3] and m show the time-dependence of the different uncertainties. The discontinuities at t = ±| 
is simply the manifestation of the fact that the functions f{t) and g{t) themselves are discontinuous 
at these times. Before the interaction was switched on, the uncertainties were having values equal 
to those for the vacuum state. During the time of nonvanishing interaction (and even after the 
interaction gets switched off!), they oscillate with frequency equal to twice that of the oscillator. 





n In 27! ijL in 



Figure 3: The time-dependences of the uncertainties A^; and Ap for mo = 1, w = 1, /o = -\/5, 



90 



V5, 



T = f and ^ = 0.01 



Figure 4: The time-dependence of A^p for the same choice of values 
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5 Implications in Quantum Optics 

In Quantum Optics a monochromatic (single-mode) coherent hght field is usually described by the 
harmonic oscillator coherent states It has also been shown that a coherent state (in particular 
the vacuum state) remains to be coherent under the FHO Hamiltonian |29]. The annihilation and 
creation operators for photons are related to the field quadratures Xi and X2 by 



a = Xi + 1X2-1 a) = Xi ~ 1X2 



(68) 



Xi and X2 being hermitian. The commutation [a, a^] = 1 translates to [Xi,X2] = |. The coherent 



i which is 



state has different uncertainties as A^i = \, A^a = \ and AxiXa = ^ Axi-A^a 
the minimum. Also, the photon count (probability for having a certain number of photons) in the 
coherent state is given by the transition probabilities of the corresponding number eigenstate and the 
profile is Poissonian. 

The FHO Hamiltonian 



Hit) 



hu{x^ + xi) + 



2h 

rriQUj' 



f{t)Xi + ^/2hmoUjg{t)X2 



(69) 



{z{t) is related with f{t) and g{t) by flTS]) ) with the effective noncommutativity between time and the 
field quadrature Xi of the form 



[t,X^ 



(70) 



will allow us to use the calculation of the previous sections. The photon count will be given by 
while the uncertainties in the field quadratures will get modified as 



Ax.(t) = 


1 + 


dmo f 
'W / 


Ax,(t) = 


1 
2 


Ah / 

J —00 


Ax,x,(t) = 


^( 

2h \ 





dr sin {2uj{t — t)) 



-fir) 



ITLqUJ 

^a'i'r) ^—fi-r) 

rriQUJ 



ujg (r) + 



-fir) 



(71) 



We further study the correlation among the photons. The time-evolved vacuum state 

\i{t 00)) = y^^Am,o\m) 



m=0 



will give 



N 



lit — 7- oo)|a^a|z(t —7- 00)) = '^^mPmfi = 1^1' 



m=l 



26 



(72) 



(73) 



N being the average number of photons in state \i{t ^ 00)). Also 

AO 



{i{t ^ oo)|aVaa|z(t ^ 00)) = ^m(m - l)P„,o = 1^1' - {lm{M'') + \i?Im{P20) • (74) 



m=2 



This, to linear order in 9, gives the 2nd order correlation among photons with zero time delay to be 
equal to (see Appendix [C]) 



^(^)(0) = 1 



S2 



ImiM'' 



Imif3se 
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Ar2 



(75) 



For the case Im{(3zi'^) < ^ fi'^^''(0) > 1, the photons try to bunch together while for Im{l3^^'^) > 
=^ fi'^^^O) < 1, they anti-bunch |24J. For the functions fjlHj) . we get 

ImiPs^ ) = go sm 2uT, (76) 

which imphes that the bunching or anti-bunching will depend only on the sign of the factor qq sin 2uT. 
For the choices taken in figures [21 [3] & IH ujT = | and hence no bunching or anti-bunching occurs. 



6 Conclusions 

In this paper we developed a formalism to compute the transitions between states of a quantum 
mechanical system with noncommutative time. We found that for a free Hamiltonian in ^ which is 
independent of time, the transitions are equal to the same for a different Hamiltonian in R^'^ found 
after the replacements (fT5|) . The time evolution of an operator and its expectation value (and hence 
also its uncertainty) can also be found in a similar manner. Specifically, for FHO the transition 
probabilities get modified and is given by ( 136|) and (jS]). The Poissonian distribution for the "vacuum 
to any state transition" also gets modified and is given by (HBl) . The study of uncertainties in position 
and momentum says that the time-evolved state is no more coherent. It gets some squeezing effect 
due to the noncommutativity, keeping the product of the uncertainties minimum. These uncertainties 
are explicitly found and is given in fl60|l . The leading order corrections in these uncertainties are 
oscillatory in time and they depend independently on the mass of the particle mo and the frequency 
of the oscillator u (note that the commutative uncertainties depend only on the product mQUj). These 
results suggest a possible modeling of the noncommutativity for the nonlinear phenomena in Quantum 
Optics. The noncommutativity results in the following nonlinear effects: 

1. The photon-count gets modified from the usual Poisson distribution. 

2. The uncertainties in the field quadratures change keeping the product minimum (the squeezing 
effect). 

3. The second order correlation function g^'^\0) gets modified producing new effects like bunching 
or anti-bunching of photons depending on the value of /m(/33^^). 

All these observations suggest that the noncommutativity produces incoherency in the otherwise 
coherent field. 

As a future work one can try to formulate the scattering process in higher dimensions and study its 
implications in quantum optics. The correspondence found in this paper between noncommutativity 
and quantum optics also encourages one to study such possibilities in other forms of time-sapce non- 
commutativity. As an example one can start with assuming the spacetime dependent noncommutative 
parameter 9 [3T1433] . 



Appendix 

A The Green's operator function 

Expanding G{t,to) as a Fourier integral we get 



1 /--t-oo 

Git, to) = Git - to) = — / du'giu')e~'-'^'-'"\ 
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where g{u}') is given by g{uj') = {fvjj' — Hq) ^. Introducing the simple harmonic eigenstate basis and 
using the completeness relation we get 

Here, ie has been introduced to avoid the pole on the contour (real axis). After finding the integral 
inside the summation using the complex analysis we get the Green's function to be 

G{t,to) = ~e{t-to)e-^'^^'-'''\ (77) 
n 

where 9(x) is the Heaviside Step Function. 



B Simplifying A^^\T2,Ti 



To simplify (l32l) what we do is to find the differential equation for A^^\t, Ti) with independent variable 
t and solve it with proper initial conditions. The differential equation has been found to be 



(^th^^ - HY^it)^ A^'\t,T,) = HrAt)A^'\t,T,), (78) 



which is a first order equation and hence A^^\t,Ti) is unique if an initial condition is given. The 
initial condition comes from the fact that A^'^\t,Ti) and A^^\t,Ti) must become identity and zero 
respectively for t = Ti, i.e., no interaction. Thus, A^^'^ {Ti,Ti) = 0. To solve the equation we define 
the Green's operator function Gmtit^to) as 

z^l- - H]l\t)\ GUtM) = Sit - to). (79) 



dt 

Now generalizing solution (1771) for the time-dependent case we get 



G.„,(t,to) = -^e(t-to)T e-^^'o^^^lS^M =-^Qit-to)A^'\t,to). (80) 



h 



I 



h 



It can be easily checked that the above expression for the Gint{t,tQ) satisfies the corresponding differ- 
ential equation. The solution for A^^\t,Ti) is then given by 

/+00 
dto GUt,to)H}Tito)A^'\to,T,), (81) 
-00 

where A^i^^^{t,Ti) is the solution of the homogeneous equation {ih-^ — HY^*{t)) A^^l^{t,Ti) = 0. The 
G-function in the expression of Gint{t,tQ) and the fact that interaction was off before Ti, with the 
initial condition for A^^\Ti,Ti) = 0, gives the initial condition for ^^^''^(t, Ti), i.e., A^^2n(Ti,Ti) = 0. 
The only solution of the homogeneous equation with this initial condition is A^i^^^{t,Ti) = 0. Thus 
we get 

/+00 
dto G,nt{t,to)H}f{to)A^''\to,T,). (82) 
-00 

Now, putting the expression of Gint(t,to) above and introducing A^^''(to,Ti) [A^^''(to,Ti)'j before 
Hffito), we get fl33|) . Here we have also used the following property of A^^\t,Ti): 

A^'\t,to)A'-'\to,T,) = A^'\t,T,). 
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C The correlation function 

An operator corresponding to the detection of a photon by a detector should be proportional to the 
annihilation operator a (say ka) [211130] • Hence if \i) is the initial state of the radiation field, the 
state after the detection of one photon is ka\i). The amplitude for going to the final state |/) is given 
by k{f\a\i). The corresponding probability is (/|a|i) p. Thus the probability of detection of one 
photon in the state \i) 

Pi = Y.\^?\m^? = |A;r5^(2|at|/)(/|a|^) = \k\^i^Wa\^. (83) 
/ / 

Similarly, probability of detection of two photons with a time delay of r is 

P2 = Y}kmHt + r)a{m\' = \k\'Y.^\a\t)a\t + r)\f){f\a{t + r)a(t)K) 

/ / 

= \k\'^{i\a\t)a)(t + r)a(t + r)a(t)|«). 

The 2nd order correlation function with a time delay r is defined as 

(2)/ N _ -P2(r) _ {i\a} {t)a\t + T)a{t + T)a{t)\i) 
^ ~ (z|at(t)a(t)K)2 • ^^^> 

For r = 

(2)fa\ {i\aKt)0'Kt)(^{t)a{t)\i) {i{t)\a^ aa\i{t)) 



(84) 



For a coherent state it can be calculated to be equal to 1. 
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